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Here we discuss some issues concerning the statistical properties of ocean surface waves. We 
show that, using the approach of weak turbulence theory, deviations from Gaussian statistics can be 
naturally included. In particular we discuss the role of bound and free modes for the determination 
of the statistical properties of the surface elevation. General formulas for skewness and kurtosis as 
a function of the spectral wave action density are here derived. 



PACS numbers: 



I. INTRODUCTION 



Recently it has been observed experimentally that 
large amplitude waves can appear on the surface of the 
ocean. There are different mechanisms that can lead to 
the formation of such events. For example, linear the- 
ory can be considered as a first candidate for explain- 
ing these waves: it can happen that for some fortuitous 
occasion, two or more waves of different lengths are in 
phase, leading to the well known constructive interfer- 
ence mechanism (linear superposition of Fourier modes). 
When the surface elevation is characterized by a large 
number of linear waves with random phases, it is possi- 
ble to estimate the probability of measuring an extreme 
event (say a wave of height larger than 2 times the sig- 
nificant wave height). This theoretical task has been ac- 
complished about 50 years ago [j| and the main result is 
that, if the surface elevation is Gaussian and the process 
is narrow-banded, then wave heights and wave crests are 
distributed according to the Rayleigh distribution (cor- 
rections due to finite spectral band- width have also been 
obtained). 

First substantial corrections to this distribution can 
be obtained if waves are considered weakly nonlinear, or 
more precisely, if for each wave component (free modes) , 
its bound contributions are included. For the narrow- 
banded case this is nothing but describing the surface 
elevation as a Stokes expansion. The more general de- 
scription of the surface elevation, valid for any spectral 
band-width, was given in a seminal paper by Longuet- 
Higgins 6]. He was able to derive the contributions to 
the surface elvation from bound waves up to second order 
in wave steepness. The numerical implementation of the 
formulas reported in the paper by Longuet-Higgins cor- 
responds to what today is called a "second order theory" 
(see 0); note that the theory includes only bound modes 
and not free modes. The presence of those Stokes-like 
contributions gives the waves the well known property of 
being positively skewed. Using the second order theory, 
[Tl| was able to include this contribution in the distri- 
bution function of the wave crests. More in particular, 
for unidirectional narrow-banded waves in infinite water 
depth, with the hypothesis that free waves are described 
by a Gaussian statistics, he derived a formula for the dis- 
tribution of wave crests (now known as the Tayfun distri- 



bution) which does bring substantial corrections to the 
Rayleigh distribution, especially if the wave steepness is 
large. It should be here stressed that the Tayfun second 
order theory predicts a Rayleigh distribution for wave 
heights (this is because second order contributions cancel 
out for wave heights). The Tayfun distribution has been 
recently compared successfully with interes ting numeri- 
cal experiments of envelope equations (see H3). Para- 
doxically, the distribution describes very well the data 
characterized by large directional spreading but under- 
estimates wave crests in the long-crested case for which 
the distribution has been derived. 

Only in the last few years it was realized that not only 
bound modes can generate deviations from a Gaussian 
statistics but also the dynamics of free waves should be 
considered in the determination of the statistical proper- 
ties of the surface elevation. More in particular, it was 
shown in and Q that the nonlinear interactions of free 
modes can substantially alter the statistical properties of 
the surface elevation. Note that in this case the statis- 
tics of free mode (without the contribution from bound 
modes) can be non Gaussian. It was also found that the 
nonlinear interactions responsible for such a deviation 
from gaussian statistics are associated with the modula- 
tional instability mechanism (also known as the Benjami- 
Feir instability) which can be thought as a quasi-resonant 
4 wave interaction. Those concepts are the bases of the 
theory developed in [j| , where a kinetic equation that in- 
cludes quasi-resonant interactions is derived. More than 
that, in the same paper, an equation that relates the 
kurtosis (fourth moment of the surface elevation) to the 
spectral wave action density is obtained. The theory in- 
cludes only the contribution to the kurtosis from free 
modes. 

The aim of the present paper is to describe a single 
theory, based on the Hamiltonian formulation of surface 
gravity waves, that can take into account both the con- 
tribution of bound and free waves. Before entering in the 
discussion we will motivate the importance of including 
the free wave dynamics in the theory by showing some 
experimental data recorded during the Marintek experi- 
ments S\ . Here we anticipate that the experimental data 
suggest that in the particular condition of long crested 
waves and large Benjamin- Feir Index, the second order 
theory is inadequate to describe the distribution of wave 



2 



crests. In the last part of the paper we will sketch the 
derivation of the formulas for the kurtosis and skewness 
which include both the contribution from bound and free 
modes. 



II. DATA FROM MARINTEK COMPARED 
WITH THE SECOND ORDER THEORY 

Here we will consider some experimental data recorded 
at Marintek in one of the largest wave tank in the world. 
The length of the tank is 270 m and its width is 10.5 
to. The conditions at the wave maker were provided by 
a Jonswap spectrum with random phases. One accepts 
this fact and then lets the phases (and amplitudes) evolve 
according to the nonlinear dynamics. Different runs were 
performed (see 8]). Here we will consider just the prob- 
ability density function of wave crests for the run char- 
acterized by a strong nonlinearity (steepness calculated 
as £ = H s k p /2 was about 0.15, with H s the significant 
wave height and k p the wave-number of the peak of the 
spectrum, computed using the linear dispersion relation 
from the peak frequency). Experimental data are com- 
pared with numerical data from a standard second order 
theory (the coupling coefficient has been taken from [l| 
and with the Tayfun distribution. 

In figure ^ we show a comparison between experimen- 
tal the wave crest distribution recorded at the first probe, 
a few wave lengths from the wave maker, and the second 
order theory. First of all, it should be mentioned that 
for the present conditions the Tayfun theory is in good 
agreement with second order theory; the experimental 
data are not so far from the mentioned distributions. 
This is not a surprise because data have been built as 
a linear superposition of Fourier modes and the gener- 
ation of bound modes is practically instantaneous with 
respect to the Benjamin-Feir space scales. The situation 
is different after the waves have evolved along the tank 
according to their nonlinear dynamics. As can be seen in 
figure |2 the Tayfun distribution (and the second order 
theory) completely underestimate the experimental data. 
This result is quite significant because it clearly tells us 
that there are some statistically stationary conditions in 
which the second order theory is completely inadequate 
to describe the probability density function of real water 
waves. 

Apparently there is no reason for the second order 
theory to fail to reproduce the experimental data; data 
(steepness and spectral shape considered) do not violate 
any of the assumption for deriving the second order the- 
ory. In the second order theory the surface elevation is 
considered as a linear superposition of non interacting 
free modes to which bound modes have been added. The 
main problem in the approach is that, as it is well known 
from Hasselmann-Zakharov theory (Q, E3)> f ree wav es 
also interact. In the next section we will consider a gen- 
eral theory which considers both the contribution from 
bound modes and free modes. 




FIG. 1: Wave crest distribution: experimental data, measured 
at 2.8 wavelengths from the wave maker, are compared with 
data from Tayfun and second order theory. 




FIG. 2: Wave crest distribution: experimental data, recorded 
at 17.1 wavelengths from the wave maker, are compared with 
data from Tayfun and second order theory. 
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III. SKEWNESS AND KURTOSIS: 
CONTRIBUTION FROM BOUND AND FREE 
WAVES 

The goal of this section is to derive formulas for the 
skewness and the kurtosis as a function of the spectral 
wave action density. The starting point for the derivation 
of the theory is the Hamiltonian description of surface 
gravity waves (see and Q ). The theory described 
below is general and can be applied to any weakly non- 
linear dispersive system. Here we will start with the fol- 
lowing Hamiltonian, truncated to four-wave interactions: 



H = I /jj aoaQdk - 



Uq\ 2 (a alal + aQaia 2 )(5 -i-2*oi2 + 



(1) 



3 / ^ 2 i,2( a o a i a 2 + aoaia 2 )<5o + i+2dfcoi2 + 



(2) * * 
K),l,2,3 a O a l a 2 a 3^0+l-2-3rffcoi23 + 



The notation is taken from |J| (for example, = a(kj), 
U^\ 2 = ^^Kkj, k l7 k 2 )). The coupling coefficients and 
their properties are reported in the just mentioned pa- 
per (the reader unfamiliar with the present Hamiltonian 
description of surface gravity waves who is willing to 
understand details of what follows is strongly advised 
to read the paper Q). Note that we have for brevity 
omitted three more integrals which contain non resonant 
four-wave interactions. Up to our desired accuracy, those 
terms do not contribute to the skewness and to the kur- 
tosis, therefore can be safely neglected. The surface el- 
evation ?7(k) and the velocity potential ip(k) are related 
to the wave action variable a(k) in the following way 
(note that the dependence on time has been omitted for 
brevity): 



77(k) 



V-(k) 



2cj(k) 



(a(k)+a*(-k)) 



2k 



(o(k)-o*(-k)). 



(2) 



(3) 



It is well known that three waves resonant interactions 
are forbidden for surface gravity waves, therefore it is 
always possible to introduce new canonical variables b 
and b* for which the Hamiltonian does not show explic- 
itly those interactions. The standard way of doing this 
consists in considering the canonical transformation, ex- 
pressed as the following infinite series: 



a = b 



A 1 ! 2&1M0- 



1-2 



dk 



12- 



J A\a b i b ^+i-2dki2- 

^CU, 2^1 ^2^0+1+2*12 



(4) 



The canonical transformation allows one to separate di- 
rectly the bound modes from the free modes: modes 6(k) 
are free modes which have their own dynamics and the 
surface elevation (which contains free and bound modes) 
can be recovered directly by using equations (0) and 10}. 
Indeed, it is easy to show that in the narrow-banded ap- 
proximation equations J2J and reduce to the Stokes 
expansion, i.e. a carrier wave plus its bound harmonics. 

Before entering into the details, here we mention that 
we will make large use of the so called quasi- gaussian 
approximation, i.e. the goal is to express the higher order 
correlators as a function of the second order correlator, < 
&162 >= N181-2, with N = N(k) the wave action density 
spectrum (< ... > indicate ensemble averages). In this 
approximation the fourth and six order correlators for 
the free waves can be writen as follows (see for example 



< 6l&2&3&4 >= ^1^2(6l-3<5 2 -4 + ^1-4^2-3) + £1,2,3,4 



< &162&3&3&I&5 > = N^NZ 

(<5l-4<52-5^3-6 + ^1-4^2-6^3-5 + 
+ <5i_ 5 (5 2 -4<53-6 + ^l-5<52-6^3-4 + 
+ (51-6^2-4(53-5 + (5i_ 6 (5 2 -553-4) 



(5) 



(6) 



-Di,2,3,4 is the cumulant, irreducible part of the fourth or- 
der correlator (for the sixth order correlator it has been 
neglected). Note that, for a gaussian process £>i,2,3,4 is 
exactly zero. Therefore, in order to describe departures 
from gaussian behavior we consider an evolution equation 
for -Di,2,3,4- This can be accomplished starting with the 
evolution equation for the free waves (the Zakharov equa- 
tion) and developing evolution equations for the second 
and fourth moment of a(k). the equation for -Di, 2 ,3.4 can 
be solved supposing that the wave action density spec- 
trum changes on a slow time scale and that at time t = 
the waves are Gaussian (the cumulant is exactly zero). 
Details can be found in [2j. The result is the following 
equation for -Di.^3.4: 



-Dl,2,3,4 =2Ti i2! 3,4<5i +2 _3_4G( Au, t) 

NxN 2 {N 3 + Ni) - (Nx + N 2 )N 3 N± 



(7) 



21,2,3,4 is the coupling coefficient for the Zakharov equa- 
tion and G(Atu,t) is given by: 



Auj 



with Acj = uj\ +UJ2 — ^3 — ^4- We now consider the statis- 
tical properties of the weakly nonlinear system described 
by the Hamiltonian in (Q. We consider the third order 
moment < ?y(x) 3 >; using the definition of the Fourier 
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transform and equation (J2J we obtain: 
< ?7(x) 3 >= J ii, 2 ,3(< aia2a 3 > + 
+ < a*a*.a 3 >)#i 2 3dfcl23 + 3 J Li,2,3(< aic^ajj > 



+ < aja^ >)5i +2 — 3^ 



L2:i 



(9) 



where 



^1,2,3 



LU1U2UJ3 



8g 3 



(10) 



Note that £1,2,3 is symmetric under the transposition of 
all subscripts. We now use the canonical transformation 
Q and insert it in © to obtain: 



< ?7(x) 3 >= J Si, 2 ,3, 4 (< M2&3&4 > + 



(11) 



< 6*626364 >)(5i+2-3-4rffcl234 



where the kernel 5*1,2,3.4 is given in [!j where also interme- 
diate steps in the derivation are reported. Here we just 
mention that 5*1,2,3.4 is a function of and 
in the canonical trnasformation. Now using equations (J5J 
and J7J) we can write the skewness as: 

< ^(x) 3 >= A J S h2X2 N 1 N 2 dk 12 + 

+ 1fi/<? T N N N 1 ~ Cos ( A ^) (12) 

+ 10/ Ol, 2,3,4^1, 2,3, 4JViJV2iV 3 

J Au 

^1+2-3-4*1234 

The final result is that the skewness can be written as 
the sum of two contributions: the first one includes only 
bound modes and the second one depends both on bound 
and free modes. 

We now consider the kurtosis. Using a similar ap- 
proach used to derive equation l|12|) , we write the kurtosis 
in the following way: 



< r?(x) 4 >= / Mi, 2 ,3, 4 (< aia 2 alal > + 



where 



+ < ala 2 a 3 a4 >)<5i+2-3-4*i234 



3 

4? 



-^1,2,3,4 = — V w l w 2^3W4 



(13) 



(14) 



Using the canonical transformation, we obtain 



< ^(x) 4 >= J 



Mi,2, 3 , 4 (< 61636*64- > 



< 61626364 >)£l+2-3-4*1234- 
#1,2,3,4,5,6(< 616263646^6* > 



(15) 



< 61626364656^ >)<5l+2+3-4-5-=fid&123456 

The kernel K 1,2, 3, 4, 5, 6, given in (9j, is a function of 
A (2) , A (3) . Note that at the same order, higher 
order terms that we have omitted in the canonical trans- 
formation should enter (those terms are reported in JjJ). 
We now apply the quasi- gaussian approximation and 
use equations JSJ and © for the fourth and sixth order 
correlator to obtain: 



< ?7(x) 4 >= 3 < r,{^Y > + 

+ 16 / Mi,2,3,4Ti 2 ,3,4iViiV 2 iV3 



1 - Cos(Aujt) 
Aw 



(16) 



<5l +2 -3-4*1234 + 12 / -^1,2,3,1,2,3^1^2^3*123 



This is the final formula for the kurtosis expressed as a 
sum of 3 contributions: the first one corresponds to a 
pure gaussian system, the second one is the contribution 
from free modes and the last one is the contribution from 
bound modes. Equations l|12|) and (|16fl are valid for long 
and short crested waves. It is interesting to note that in 
the limit of long-crested waves and in the narrow-banded 
approximation the integral including free modes becomes 
proportional to the Benjamin-Feir Index (see [3j for de- 
tails), while the integral including bound modes become 
proportional to the steepness squared. 

As a conclusion we may state that the main contribu- 
tion of the paper is the derivation of equations (|12fl and 
(|16fl . Those equations, obtained using the Hamiltonian 
formulation of surface gravity waves, can be viewed as a 
generalization and unification of the theory developed by 
Longuet-Higgins in 1963 on bound modes and the theory 
developed by Janssen in 2003 on free modes. 
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